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Abstract 

We calculate the inclusive small-x valence quark production cross section in proton- 
nucleus collisions at high energies. The calculation is performed in the framework of the 
Color Glass Condensate formalism. We consider both the case when the valence quark 
originates inside the nucleus and the case when it originates inside the proton. We first 
calculate the cross section in the quasi-classical approximation resumming the multiple 
rescatterings with the nucleus. Then we include the the effects of double logarithmic 
reggeon evolution and leading logarithmic gluon evolution in the obtained cross section. 
The calculated nuclear modification factor for the stopped baryons exhibits Cronin en- 
hancement in the quasi-classical approximation and suppression at high energies / rapidities 
when quantum evolution corrections are included, providing a new observable which can 
be used to test Color Glass physics. 
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1 Introduction 



One of the most striking features of experimental data in nuclear collisions at high energies 
is the large stopping suffered by the nucleons entering the collision [H 13 El- This stopping 
translates into a net baryon number transfer from the forward to the central rapidity region, 
and has already been observed in proton-nucleus and nucleus-nucleus collisions at CERN jH |2] 
and in Au + Au collisions at RHIC 0. The description of this phenomenon has been previously 
approached in the framework of non-perturbative models for hadron collisions El El 13 El • In 
in El El m El the baryon number transfer along a large rapidity gap is attributed to topological 
excitations of the gluon field. 

In this work, following IH], we consider the problem of baryon stopping from a purely 
perturbative perspective in which the carriers of baryon number are the energetic valence quarks 
entering the collision, that loose much of their longitudinal momenta in the collision and are 
driven towards the central rapidity region via hard gluon emissions. Valence quark production 
in the scattering of a proton on a nucleus is calculated in the framework of the Color Glass 
Condensate (CGC) formalism |iniinill2llIH|llllliniliniinilIH|. The effects of strong gluonic 
fields in the small-x hadronic or nuclear wave functions on baryon number transfer were first 
considered by one of the authors with collaborators in The valence quark distribution 
at small Bjorken x was first studied in |U] in the quasi-classical approximation of McLerran- 
Venugopalan model JT] . Quantum evolution corrections come into the obtained valence quark 
distribution through a non-linear evolution equation which was derived in The linear part 
of the evolution equation from [9j corresponds to the well-known double logarithmic reggeon 
evolution equation from EDI 1^ 1221- It is an interesting property of quark evolution that 
at the leading order each power of the coupling constant is accompanied by two powers of 
the logarithm of center-of-mass energy Ins, or, equivalently, of the logarithm Inl/x of inverse 
Bjorken x ^21 EOl 1^ I22|- The resummation parameter is thus In^l/x, unlike, say, the 
BFKL equation [2n! , which resums powers of as In 1/ x. Due to that property one may expect 
the small- a; evolution corrections to manifest themselves earlier (i.e., at larger values of x 
and smaller rapidities) in observables related to valence quarks. The non-linear part of the 
evolution equation derived in P is due to gluon evolution effects, similar to those resummed in 
[13 Cni im CHI • While the linear reggeon evolution [121120112111221 tends to increase the small- a; 
valence quark distribution, the effect of nonlinear evolution equation [9 is to reduce the valence 
quark distribution at very small-x as compared to the linear evolution predictions. 

Here we continue the work done in ^ by calculating the inclusive production cross section for 
small-a; valence quarks in pA collisions. We perform the calculations both in the quasi-classical 
approximation [TTl [T21 and including quantum evolution from P . We use the techniques for 
calculation of inclusive cross sections developed in [^EH] (for a review see |2E1)- The calculation 
of inclusive gluon production in pA collisions was performed in [^ EZj in the quasi-classical 
approximation and in [23 the effects of quantum evolution were included in the cross section. 
Our calculations will proceed along similar lines. In Section 2 we derive the expression for 
valence quark production in the quasi-classical approximation resumming multiple rescatterings 
of the incoming proton and produced valence quark on the nucleons in the target nucleus. Here 
one has to separately consider two cases: the valence quark may originate either in the proton 
or in the nucleus. The stopping of valence quarks in these two cases is clearly distinguishable 
experimentally: the rapidity distribution of stopped valence quarks from the proton scales 
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as ~ e^~^ while the valence quarks from the nucleus are distributed as ~ e~^. (Here the 
nucleus has rapidity 0, the proton has rapidity Y and the quark is produced at rapidity y.) 
We calculate both contributions in Sect 2.1 and study their properties in Sect 2.2. Our results 
here are complimentary to the valence quark production calculated in [28\: the authors of [2H] 
considered production of a hard valence quark which experiences no recoil and is produced in 
the fragmentation region, while here we are interested in production of a soft valence quark far 
away (in rapidity) from the fragmentation region. 

As was shown in [SHI IHOl ISH IHS IHHl IHH tHEl IHEl ini| for gluons, multiple rescatterings in 
the nucleus lead to enhancement of high-py particle production which is similar to the Cronin 
effect observed in jSFj. Analyzing the nuclear modification factor resulting from the derived 
quasi-classical cross section for valence quark production, we observe that valence quarks also 
exhibit Cronin enhancement, as demonstrated in Fig. 

We continue in Section 3 by including the effects of non-linear reggeon evolution from 12] 
into the valence quark production cross section. The cases of valence quark originating in the 
proton and in the nucleus are considered separately in Sections 3.1 and 3.3. Similar to the case 
of gluon production [23], the evolution between the produced quark and the projectile turns 
out to be linear and is given by the reggeon evolution equation ^H] if the quark comes from 
the proton and by the BFKL equation [23j if the quark comes from the nucleus. The evolution 
between the produced valence quark and the target nucleus is always non-linear. 

As was originally observed in 1211 110] for gluon production, small-x evolution introduces 
suppression in the nuclear modification factor at high energies/rapidities, turning Cronin en- 
hancement into suppression. These predictions were confirmed by the forward rapidity particle 
production data produced by RHIC d + Au experiments jiOl HH 1121 UHl HI]- In Section 3.2 we 
study the effects of non-linear reggeon evolution on valence quark nuclear modification factor 
for the case when the valence quarks are produced from the proton. We observe suppression 
of nuclear modification factor which increases with rapidity /energy similar to the case of the 
gluons, as shown in Fig. [71 Therefore, if baryon stopping is dominated by perturbative mech- 
anisms, we expect the nuclear modification factor for net baryons to be suppressed at forward 
rapidities in d+Au collisions at RHIC. Thus the nuclear modification factor for stopped baryons 
may serve as a new test of Color Glass Condensate (CGC) dynamics in nuclear collisions, or, 
at least it may help determine the relative contributions of perturbative and non-perturbative 
baryon stopping. 

We conclude in Section 4 by discussing the applicability of obtained results and the prospects 
for experimental testing of the derived formulas. In order to compare our results to the actual 
experimental data produced at RHIC and to make predictions for the LHC, as was done in 
|45| [3S1 IS [H]; one should implement the fragmentation mechanism which converts valence 
quarks into hadrons. This goes beyond the scope of this paper and is left for future work. 
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2 Valence Quark Production in the Quasi-Classical Ap- 
proximation 

2.1 Derivation of the Expression for the Production Cross Section 

We are interested in calculating the soft valence quark spectrum in proton-nucleus collisions 
at high energies. In what follows we will assume that the proton is moving ultrarelativistically 
in the light-cone 'plus' direction and, consequently, has a large Pi momentum component, 
whereas a nucleon in the nucleus is moving in the light-cone 'minus' direction with a large P2 
momentum component. We employ Sudakov parameterization for the light-cone components 
of the momentum of the the produced valence quark: k = {api , l3p2 ,k). In the center of mass 
frame of the collision, the relation of the dimensionless Sudakov parameters a and j3 to the 
rapidity of the produced quark y and the total rapidity of the collision Y is 

« = ^e^->'/2, P=\^e-y^''/\ (1) 

where s = 2p^ is the square of the center-of-mass energy of the collision and the proton has 
rapidity y = Y while the target nucleus has y = 0. The calculation below will be performed 
in the light-cone gauge of the proton, A'^ = 0, which for the nucleus moving in the 'minus' 
direction is equivalent to the covariant gauge, d ■ A = 0. This choice completely determines 
the set of diagrams that are relevant for our calculation as well as their physical interpretation: 
As it was discussed in detail in [211 EI], in this gauge the classical current associated with 
the nucleus, J~ , remains unchanged throughout the collision and the nuclear effects become 
manifest in the form of multiple scatterings between the nucleons in the nucleus and the rest 
of the system. 

As usual at high energies, the resummation of the multiple rescatterings is performed in 
the eikonal approximation, in which the energetic projectiles (quark or gluon) remain at fixed 
transverse coordinate as they propagate through nuclear matter. Hence the calculation is done 
in coordinate space. Besides, we restrict the interactions with the nucleons in the nucleus to 
two-gluon exchanges, in the spirit of the quasi-classical approximation. More technically this 
implies that the resummation parameter is a'^A^^^, with A the atomic number of the nucleus. 
We are assuming that the energy of the collision is large enough for an eikonal description of 
the multiple rescatterings to be appropriate, but not large enough for quantum corrections to 
become important. Under these two approximations (eikonal limit and quasi-classical descrip- 
tion) the effects of multiple rescatterings, and therefore all the nuclear effects, can be recast in 
the form of Glauber-Mueller propagators PU ES] ■ 

We divide the cross section for soft quark production into two terms: the first contribution, 
labeled (a), corresponds to produced valence quark coming from the proton, while the second 
one, labeled (b), corresponds to valence quark coming from one of the nucleons in the nucleus 
via interaction with a soft gluon in the proton wavefunction (see below): 

da da"- da^ , , 

+ -HTi:.- (2) 



(P'k dy d'^k dy d^k dy 

Importantly, the formation time associated with the soft valence quark production mecha- 
nism is, in both cases, much longer than the typical time scale for the nucleus-proton interac- 
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tions, which can be considered instantaneous. Moreover, the diagrams in which the soft valence 
quark emission happens during the interaction with the nucleus are suppressed by powers of 
the total energy of the collision and, therefore, can be neglected in the high energy limit. This 
allows us to deal with time-ordered diagrams or, more precisely, the calculation is more con- 
venient (and will be done) in the framework of time-ordered light-cone perturbation theory 



nucleus 



proton 



nucleus 




i>n> 




B 



nucleus 



Figure 1: Diagrams corresponding to the first contribution to the cross section as considered 
in the text, where the produced valence quark originates in the proton. 

First we are going to find the contribution to the cross section in which the produced soft 
valence quark originates in the proton wavefunction. The contributing diagrams are shown in 
Fig. ^ The discussion here will proceed along the lines of the calculation done in of the 
gluon production in pA collisions. The physical picture is the following: one of the fast valence 
quarks in the proton has to emit a hard gluon in order to loose most of its large longitudinal 
momentum and, therefore, become soft. This splitting is given by the soft-quark wavefunction 
of the proton at the lowest order in perturbation theory (see appendix A): 

r.A^,z, a) = ge[l - a A] ^ , (3) 

where z is the transverse coordinate of the hard gluon (which is equal to the transverse coordi- 
nate of the fast quark), x is the transverse coordinate of the soft quark, a is the helicity of the 
fast quark, a and A are the color and polarization of the hard gluon. The helicity of the soft 
quark is equal to the one of the fast one due to helicity conservation for massless quarks. 

Analogously to what was done in [^j, it is convenient to introduce the soft quark distribution 
of a proton. In order to do so we multiply the above wavefunction in Eq. Q by its complex 
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conjugate evaluated at transverse position y, average over the helicity of the initial quark and 
sum over the polarization and color of the hard gluon, getting 



1 a 



as aNc {x-z)-{y- z) 



(4) 



d'^z dy 



Air 2 



27T {x — zy{y — z_) 



2 ' 



where dg = asCp/T^- The factor of Nc on the right hand side of Eq. (j3]) accounts for the three 
valence quarks in a proton. All the dynamics corresponding to the soft-quark emission in the 
diagrams we are going to calculate is included in the distribution in Eq. (jH). 

Our next step is to resum the multiple rescatterings of the quark-gluon system in the 
nucleus. The diagram shown in Fig. corresponds to the case in which the incoming proton 
has already developed a soft valence quark component in its wavefunction before interacting 
with the nucleus both in the amplitude and in its complex conjugate. As is explicitly shown in 
Fig. in this case only the interactions between the nucleons and the soft-quark line survive. 
The interactions between the nucleons and the hard gluon vanish due to real- virtual cancellation 
[21]. This can be understood in the following terms: since we want to keep the momentum of 
the produced soft valence quark fixed, in the transverse coordinate space that we are going to 
perform our calculation in, its transverse coordinate in the amplitude x_ has to be different from 
the one in the complex conjugate amplitude y j^EE]- The transverse momenta of all other 
particles in the final state are integrated over (inclusive production) making their coordinates 
in the amplitude and in the complex conjugate amplitude equal to each other. Therefore, one 
can freely move the t-channel gluon lines connecting the nucleons with the hard gluon across 
the cut without changing either the color factor or the momentum of the produced quark, but 
picking up a relative minus sign in the way, which provides the cancellation. 

To calculate the surviving contribution one can proceed either by direct evaluation of the 
diagrams or, alternatively, by making use of the crossing symmetry property |5H 1^ in the 
following way: mirror-reflecting the soft quark line in the complex conjugate amplitude of 
Fig. n\ with respect to the final state cut one is left with the diagram corresponding to the 
amplitude of a quark- antiquark dipole of transverse size |x — ?/| multiply scattering off a nucleus. 
The amplitude for this process is jSS] 



is the saturation scale of a quark dipole, p is the nuclear density (normalized to A) and T4(&) is 
the nuclear thickness at impact parameter h in its own rest frame, which for a spherical nucleus 
of radius R is equal to T^(&) = 2a/ W' — 6^. 

The diagram in which the soft quark is emitted before the interaction with the nucleus in 
the amplitude, but after the interaction in the complex conjugate one is shown in Fig. 03. 
The interactions between the nucleus and the hard gluon line in the amplitude do not cancel 
in this case since, to move the gluon line attached to the fast quark in the complex conjugate 
amplitude across the cut, one is forced to 'jump' over the emission vertex, which would modify 
the momentum of the produced quark, ruining the cancellation. To explicitly include both 
the interactions with the s-channel quark and gluon lines, we do not connect the t-channel 



N{x,y) = 1 - e"l-"^l^'^»9^°(^/l-"^l'^)/'*, 



(5) 



where 




(6) 
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Figure 2: Mirror reflection of tlie diagram in Fig. ti(2) is tlie emission time in tlie amplitude 
(complex conjugate amplitude). The interaction time is r = 0, represented in the figure by a 
vertical line and the final state cut is denoted by r = +oo. 

gluon lines to any particular line in Fig. CJ3. Again, mirror-reflecting the complex conjugate 
amplitude with respect to the cut generates the configuration represented in Fig. |21 There, 
the fast valence quark line in the complex conjugate amplitude turns into an antiquark line 
overlapping the hard gluon line in the amplitude at transverse coordinate 0. Since the whole 
system is color neutral, the zero-size antiquark-gluon system has to be in the 3 representation 
of SU{3), which together with the quark line located at the x coordinate is equivalent to a 
quark dipole multiply rescattering off a nucleus. Therefore, the factor entering the production 
cross-section due to interactions in Fig. ^3 is 

N{x, z) = l- e-l^-^l'«'^ ln{l/k-£|A)/4_ 

The contribution of the diagram complex conjugate to Fig. [Tp (not shown) is calculated in 
an analogous way yielding 

N{y, z) = l- e'^y-^^'^'" Hi/\y~m/\ (8) 



The diagram in Fig. gives a zero contribution to the production cross-section, since, due 
to real- virtual cancellations, all the interactions with the target cancel, leading to = 0. 

Finally, adding up the contributions to the production cross sections of the diagrams in 
Fig. HI gives 



da" 
(Pk dy 



(2n) 



d'^x d'^y d^z e 



ik{x—y) 



dn'^ 



X 



-ix-yfQ-i^lnil/\x~y\A)/4 _ ^^(x~z)^Q%lnil/\x-z\A)/A _ ^-(y-^fQ%Hl/\y~z\A)/i ^ ^ 



d'^z dy 



(9) 



Our next step is to calculate the diagram in which the produced valence quark originates 
in the nucleus wavefunction, an example of which is shown in Fig. |21 Before starting this 
calculation, let us discuss the diagrams that do not contribute to the cross-section in this 
case. First, the emission of the hard gluon cannot happen via interaction with any other 
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nucleon in the nucleus. In its rest frame the nucleus is a dilute system, i.e, its constituents 
nucleons are spatially separated from each other. In a frame in which the nucleus is fast 
moving in the light-cone 'minus' direction, diluteness of the nucleus translates into the fact 
that the nucleons are strongly localized around fixed well differentiated coordinates: < 



< x~^. Moreover, the gluon field of the nucleus as given by the classical equations of 



motion in the = gauge is equivalent to the gluon field in covariant gauge and is strongly 
localized around the x"*" coordinates of the nucleons jHSlE]: A ~ ^iS{x~^ ~^t)- Therefore, the 
interactions between the nucleons are suppressed by powers of A. In the leading order in a'^A^^^ 
calculation, corresponding to the quasi-classical limit used here, the nucleons can be assumed 
as not interacting with each other. On the other hand, those diagrams in which a nucleon 
exchanges one gluon with one of the quarks in the proton and subsequently emits a hard gluon 
are suppressed by a power of center-of-mass energy at high energies in A'^ = light-cone gauge 
and are also subleading. This leaves us with a single possibility for the nucleon's valence quark 
to be produced at central rapidity: it may only happen via interaction with a soft gluon in the 
proton's wavefunction, as shown in Fig. El 



proton 



1> 




X 




nucleus 



Figure 3: Valence quark production diagram corresponding to the second case considered in 
the text, where the produced soft valence quark originates from the nuclear wave function. 

The physics is clear: the proton has to emit a soft gluon before its interaction with the 
nucleus. This emission is given by the soft-gluon wavefunction of a valence quark 

<.fel) = 29*»^<p^. (10) 

The emitted soft gluon undergoes multiple rescatterings until being absorbed by one of the 
valence quarks of the zth nucleon at the longitudinal coordinate u. This way, the valence quark 
is kicked out from its original light-cone trajectory towards the central rapidity region and 
is likely to rescatter on the remaining nucleons in the nucleus. A space-time picture of the 
collision is depicted in Fig. |31A.. Integrating over the longitudinal coordinate u and summing 
over all nucleons in the nucleus one gets the following factor: 

Jo 
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The first exponential term in Eq. accounts for tlie propagation of a gluon dipole until the 
longitudinal coordinate u (see Fig. 0f3) with 

Q% = 4nalpTM (12) 




X-oo 



proton + nucleus 

^1 



Figure 4: Space-time picture of the diagram in Fig. |21 (left), and its mirror reflection (right): 



the gluon dipole saturation scale. The factor V{x, y) is the Fourier transform of the Gq — > qG 
cross-section in the high energy limit (see Appendix B for a detailed derivation): 



V{x, y) = N, I dH S-^-y-^ ^^^^^ ^ na^^Gpl In ■ ^ 



dH 



s {x — y^A? ' 



(13) 



where s = 2k^p2 = \k\^/s e'^~~ . The factor of entering the definition of V accounts for the 
three valence quarks in a nucleon. Finally, the second exponential term in Eq. (jllj) corresponds 
to the propagation of a quark dipole from u till the edge of the nucleus L, which is equal to 
Ta(6) evaluated at a fixed value of the impact parameter. Thus, the contribution to the cross 
section of the diagram in Fig. El is: 



da'' 



d^kdy (27r)2 



d^xd'^yd^ze'^^^-y) 



dn9 
d^z dy 



where 



X 



p du e~^-~^^^'^'^3 ln(l/|£-?/|A)n/4L y(^^^ ^-{x-y)^Q% \nil/\x-y\A){L-u) /4L 

.Jo ~ 

dn^ as {x-z)-{y- z) 



d^zdy TT {x — zy{y — zy^ 
is the soft gluon distribution of the valence quark. 



(14) 
(15) 
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Rewriting a in Eq. Q and /9 in Eq. ()14|) using Eqs. (0), performing the integral over the 
longitudinal coordinate u in Eq. (jl4j) and shifting the integration variables in the following way 



X — z 



y- z = y, 
z = b, 



(16) 
(17) 
(18) 



the final result for the soft valence quark production cross section in the quasi-classical approx- 
imation can be rewritten as: 



da 



da'' da'' 
+ 



d^kdy d^kdy d'^kdy' 



(19) 



with 



da"" _ 1 a,Nc y_Y/2 \M_ f ^2 .2 ik(x-y) ^ • y 
d^kdy ~ (27r)2 27r ^Jdxdyd be-- - 



X 



-{3I-?/)'Q^,ln(l/|ir-j/|A)/4 



-ir'QLln(l/|2;|A)/4 



y'QL ln(l/|y|A)/4 



+ 1 



(20) 



and 



da"" 



1 2a, C 



F e 



d^kdy (27r)2 vr 1 - ^ ^\k 



d^xd^yd%e'^^^-y) 



x-y 



x^y"^ {x — yy 



X 



-ix-y)^Q-i^ln{l/\x-y\A)/4 _ f,-{x-y)^Q%Hl/\x-y\A)/4 



(21) 



In arriving at Eq. ^IT\ we have used the definition of the saturation scale from Eqs. © and 
jni) with Ta(6) replaced by L. 

Equations (fT^ . (plj) and (PT|) provide the result for the small-x valence quark production 
in pA collisions at high energies in the quasi-classical approximation. 

2.2 Properties of the Quasi-Classical Cross Section 

In this Section we explore some of the properties of the valence quark production cross section 
derived in the previous Section. 

First let us consider the case of kx <^ Qs, where fcy = |^|. For momenta of the produced 
quark much smaller than the saturation scale, fcy <^ Qs, the dynamics is governed by the 
strong nuclear effects. In this region the integrals in Eqs. (j^UI) and ()21|) are dominated by large 
distances and setting the logarithms in the exponentials to a constant value (more precisely to 
1) is a good approximation. This way the integrals are analytically doable, yielding 



da'' _ d,N, e^-^/^ ^ 



d^kdy 27r ^/s 



1 2 uijci'^ 
rrr + Trre - '^"^ 



sq 



Ei 



+ In 



Q 



4 

sq 



4rA2 



(22) 
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and 



-y+Y/2 ^ 



d?k dy 



d^b { In 



r 



0% 

4A2 



sq 



+ lE 



+ In 



r 



'g2 



+ In 



,2 H 



(23) 



where Ei is the exponential integral function and F is the incomplete gamma function. 

On the other hand, for very large transverse momenta, fcy 3> Qs, the nuclear effects are 
absent (tend to disappear) and the perturbative result should be recovered. Expanding the 
exponentials in Eq. ^ and Eq. (fT^ one gets the following expansion in powers of Qsqlkx 
("twists"): 



da'' 
d^k dy 



2ti 

+2 In 



4A2 



2j Qsq 



d'h 



U3 



In 



4A2 



+ 27£; - 1 



+ 



Q 



2^" 4A^ 



(37i5 - 4) + 67I - + 



29 



(24) 



and 



da^ 



d'^k dy 



271 



d% 



Ql 



U3 

IXirp 



k^ 

In ^ + 27s 



4A2 



+ 



Q 



3, A;\ 



l) + 3(27 



47i? + 1) 



+ ... 



+ 1 



3 



8 4A2 



(25) 



Remarkably, at large transverse momentum and fixed rapidity the valence quark production 
cross section scales like ~ which is to be compared with the ~ l//c|. scaling of the inclusive 

gluon production cross section. As is well-known, the quark production cross section turns out 
to be more sensitive to the ultraviolet. 

To study the transition between the small and large transverse momentum regions and to 
have a quantitative measure of the nuclear effects we introduce the nuclear modification factor, 
equal to the ratio of the valence quark production cross section in pA collisions over the one in 
pp collisions scaled by the number of binary collisions 



RpA{y, kr) 



cPkdy 



A 



(26) 



d^kdy 



The cross section for proton-proton collisions can be obtained as a particular case of our cal- 
culation, just keeping the first terms in the twist expansions of Eqs. fl^ and and setting 
A = 1. We get 



A 



daPP asNc cosh(?/ - Y/2) 



d^k dy 



vr 



Qi 



4A2 



(27) 



where we have used Q^^ ~ A^^^ and the fact that, for a cylindrical nucleus, the integral over 



impact parameter gives just a factor of the transverse area of the nucleus, Sa ~ A^/^. In what 
follows we will restrict our discussion about RpA to the central rapidity region, y = Y/2. The 



10 



< 2 




J I I I I I I I I I I I I I I I I I I I \ I I I I L. 

2 4 6 8 10 12 14 

(GeV) 



Figure 5: Nuclear modification factor at central rapidity for Q^^ = 1 GeV^ (solid) and Q^^ = 
2 GeV^ (dashed). The cutoff is A = 0.1 Q^g. 

limiting behavior of RpA at momenta kx ^ Qsg can be obtained by expanding Eq. ()22|) and 
Eq. in powers of kx/Qsg, yielding 

1.2 7,4 / Q2 nr AT \ 

R^y ^ Y/2, « Q„) . ^ + ^ + ,„ ^ + ,n ^ - 2) + . . . < 1, (28, 

indicating that deep in the saturation region the valence quark production is suppressed due 
to nuclear effects and that RpA increases with increasing momentum. On the other hand, from 
Eq. ()24j). Eq. and Eq. ^I7\ we derive the limiting behavior at large momentum, ^ Qsg, 

RpAiy = Y/2, kr » Qsg) ^ 1 + y ^ In^ A_ + . . . > 1, (29) 

which shows that RpA approaches unity from above at asymptotically large transverse momen- 
tum. Therefore RpA must reach a maximum value at intermediate momentum kx ~ Qsg, as it 
is shown in Fig.El where we plot the nuclear modification factor RpA at mid-rapidity {y = Y/2) 
built up by evaluating the integrals in Eq. ()20|) and Eq. ()2H) numerically for two different values 
of the quark saturation scale. 

We therefore conclude that, in the framework of the quasi-classical approximation, the nu- 
clear modification factor for small-x valence quark production in the central rapidity region is 
less than 1 at kx < Qsg and has a Cronin enhancement at high kx ~ Qsg- An equivalent (anal- 
ogous) conclusion could be drawn for values of rapidity outside the central region. However, 
at more forward rapidities the nuclear wavefunction is probed at smaller values of Bjorken-x, 
where the effects of quantum corrections become important and may significantly alter this 
conclusion. This is known to be the case for inclusive particle production in pA collisions. 
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for which the Cronin enhancement observed at central rapidity is completely washed out by 
quantum evolution, turning it into a relative suppression at forward rapidities, as was predicted 
for gluon production in [221 EOI and experimentally confirmed in j56j . 



3 Including Quantum Evolution 

In this Section we are going to introduce the small-x evolution corrections into Eqs. ()2()|1 and 
(1211) . Since we are interested in valence quark production, the small-x evolution corrections to 
Eqs. (j2ni) and ()21|) will include both the contributions from non-linear evolution equations for 
gluon evolution CH US CHI El dl and for reggeon evolution [131201121112213 (with quarks 
in t-channel). 



3.1 Valence Quarks from the Proton 

We begin by considering case (a) above, in which the valence quark originates in the incoming 
proton. For the quasi-classical case the diagrams of this process are shown in Fig. ^ with the 
expression for the production cross section given by Eq. ()20|). We begin by rewriting Eq. ()2U|] 
with the help of Eq. ^ as 



d?kdy (27r)2 27r j ^ \x-h\'^\y-W 

X [N{x,h) + N{y,h)-N{x,y)\. (30) 



Inclusion of quantum evolution corrections will proceed along the lines of [25, (see also 
IHZl EH Eniinni ) ■ The evolution is divided into two components: emissions of gluons harder than 
the produced valence quark (with rapidities greater than y) and the emission of gluons softer 
than the valence quark (with rapidities smaller than y). The emission of harder gluons leads to 
diagrams with quarks in the t-channel, leading to reggeon evolution equations resumming double 
leading logarithms of energy (powers of a, In^ s) [IJl |2nil2Ill221in|- Here we will consider the 
reggeon evolution in the transverse coordinate-space formalism developed in [H]. The building 
block of that evolution is a decay of a valence quark into a hard gluon and a much softer valence 
quark. 

Similar to [53^ one can show that the incoming valence quark in the proton may evolve by 
emitting hard gluons only before the interaction with the target (at light cone times r < 0). 
Only one such emission in the final state after the interaction (at r > 0) is allowed: more final 
state emission would not generate longitudinal logarithms needed for the double- logarithmic 
approximation to the reggeon evolution employed here and in Pj. Thus the effect of harder 
gluon emissions is to "prepare" a valence quark which then will serve as an incoming quark in 
Fig. ^ and, in turn, may split into a soft quark and a hard gluon either before or after the 
interaction, as shown in Fig. ^ 

To include the effects of harder gluon emissions we first have to consider a slightly more 
sophisticated model of a proton. Up until now we have modeled the projectile proton by an 
incoming valence quark. To be able to construct a transverse coordinate analogy of the reggeon 
evolution equation, as was done in [H], we will model the incoming proton by a dipole consisting 
of a quark at a transverse space position z_i and an anti-quark at Zq- The dipole size z^q = z^ — Zq 
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is comparable to the proton's transverse size ~ 1/Kqcd- If the quark in the dipole carries a 
fraction ai of its longitudinal (plus) momentum component, the total rapidity interval in the 
dipole-nucleus scattering is defined by F = ln(s zIq) and the rapidity of the valence quark in the 
wave function is yi = ln(ai s zfo) [H], where s is the center of mass energy of the proton- nucleus 
system. Eq. (jHUj) now becomes 



{ziQ,ai) = — d xd yd zie 



d^kdy ' (27r)2 27r ai j l^l — ^ip |y — ^ip 

X [N{x,z^) + N{y,z,) - N{x,y)\ , (31) 

where a is defined by Eq. (J]). Indeed to obtain the valence quark production cross section 
one has to convolute the cross section in Eq. ()31|) with the dipole's light cone wave function 
squared integrating over z^^ and ai. Such wave function is well-known for the case of deep 
inelastic scattering (DIS) on the nucleus jHI]. Thus, rigorously speaking, our discussion in this 
Section will be most relevant for DIS. However, one could use the formula we obtain below for 
proton-nucleus scattering by putting 210 ~ 1/Kqcd in Y and by using ai ~ 1/3. 

We denote by r{zj^,z^,ai;z,a)/ai the probability density of finding a soft valence quark 
at transverse coordinate z carrying longitudinal momentum fraction greater than or equal to 
a in the wave function of a dipole with the initial longitudinal momentum fraction carried 
by the quark being ai. The quantity r obeys the following evolution equation in the large- A'^c 
limit 

as /-"iminii.^o^/^lj ^2^^ 
r{,ZjQ,z^,ai;z,a) = 5 {z- z^) + / ^ ?^(^i, ^2> ^> (32) 

Eq. ()32|1 is illustrated in Fig. IHl One can immediately see that Eq. (jH^ is a linearized version 







Figure 6: Graphical representation of the evolution equation (j!?^ in the text. 

of the nonlinear Eq. (43) in [Oj. In Fig. IHlthe incoming quark 1 in the dipole splits into a hard 
gluon 1 (denoted by a double line on the right hand side of Fig. ^ and into a soft quark 2. 
The subsequent evolution continues in the new dipole 12. In the non-linear evolution equation 
derived in ^ the evolution could also continue in the original dipole 01 in the sense of gluon 
(BK) evolution [THIIIHI. However, for an inclusive quantity, such as valence quark production 
considered here, all further interactions of dipole 01 cancel by real- virtual cancellations. This 
is similar to the case of inclusive gluon production cross section considered in [2^1 • 

The upper limit of the integral over a' (the longitudinal momentum fraction of the soft 
valence quark at ^2) in Eq. (jH^ results from two considerations: on the one hand, according 
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to the rules of light-cone perturbation theory [221 while, on the other hand, rapidity 

ordering of the produced dipoles requires 



z 



a' < ai 4^ (33) 



z 



as the rapidity of the soft valence quark in dipole 12 is y' = ln{a'szf2)- (To understand this 
last formula for y' remember that, in transverse momentum space, y — {Y/2) = lr^{a^/s/kT) 
and Y = \n{s/k^) leading to y = \n{as/k^).) 

For the dipole-nucleus scattering the effect of harder gluons can be included in the produc- 
tion cross section by 

UiO'"i) = 77^— ^— ! d^x(fy(fzi£ze'-'^-~y-U\zQ,z^,ai]z,a)^~ ~^ ~^ 
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(Phdy^-'""' " (27r)2 2n ai J " ' ' '\x-z\^\y-z 

X [N{x, z) + N{y, z) - N{x, y)] , (34) 

where we use y = ln(Q; s/A;|,) to relate a and y. 

Inclusion of the evolution effects due to emissions of softer gluon having rapidity less than 
y is quite straightforward and goes exactly parallel to (23 IIHI ISH IHHl IHEl lEHj • To include these 
softer emissions we need to replace 

N{x,y) ^ N{x,y,y) (35) 
in Eq. ()34p . where now N{x, y, y) has to be found from the non-linear evolution equation |15t lTE] 



dN{x^,x^,Y) ^ a^Nc [ 2 ^01 
dY 27r2 / a;2 



d^X2 ^2 2 [^(^0>^2>^) +^(^2,^l>^) ~N{Xf^,X-yX) 
20 ^21 



with the initial condition 



N{x^.X2,Y)N{x2,x,,Y)] (36) 



A^(a^o,aii,F = 0) = iV(xo,Xi). (37) 



The final answer for the valence quark production cross section including quantum evolution 
with the valence quark coming from the projectile proton reads 



^^^(^lo^tti) = (^~y2~|~~ j d^xd'^yd'^zid^ze'-^- r{zQ, z^,ai] z,a) _ ^ ^ 

X [N{x, z, y) + N{y, z, y) - N{x, y, y)] . (38) 

Once again we remind the reader that here Y ^ yi = ln{aiS Zq^) and y = ln{a s / k'^) . When 
integrating over z in Eq. ()38|) one has to keep in mind that rapidity ordering requires that the 
rapidity of the quark at transverse coordinate z, given hj y = ln[a s (^ — ^o)^]' limited to the 
interval y < y < Y , which imposes a constraint on the range of ^-integration. 

The essential ingredients of Eq. (jHH|l are the Reggeon exchange amplitude r obeying the 
evolution equation (jH^ which sums up leading double logarithms ag In^ s, and the amplitude 
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obeying Eq. ()36p summing up usual leading logarithms Ins. We remind the reader that 
the fact that we are using the amplitude r obeying a linear evolution equation (j32j) to describe 
the evolution between the produced valence quark and the projectile is not an approximation 
but an exact result of cancellations of non-linear Reggeon evolution corrections from jU] in that 
rapidity region. 

If one is interested only in rapidity dependence of the cross section, integration of Eq. 
over kx yields 



da"" , asNc a f ^ ,2 ,2 f ^ 1 

2.10! ctij = I a X a Zid zrizr., z-., ai, z, a) - 

TT ai J \x — z\ 



dy 



N{x,z,y), (39) 



where we have also made use of the fact that N{x,x, y) = 0. In the leading double logarithmic 
approximation we find a from y ^ ln{a s/Q'^{y)) where Ql{y) is the quark saturation scale of 
the nucleus at rapidity ?/, which we use as the typical transverse momentum of the produced 
quarks. Ql{y) is rapidity-dependent due to the effects of nonlinear evolution equation (jHS)). 



3.2 Properties of the Cross Section with Evolution 



First of all, let us point out that Eq. (j32j) can be solved analytically if one averages it over the 
directions of Zq^ (or Zgi)- Defining 

p = z-zo (40) 



we can rewrite Eq. (|32|) as 

r(zoi,ai;p,a) = 5 (p-^10) + — / — ^(^21, a ; a), (41) 

2t\ a zi2 

which, after angular averaging, becomes 

r(zoi,ai;p,a) = 7^ b{p - z^q) + — / — r(2;2i, a ; r, a). (42) 

2 7rp 2t[ a' ZI2 



Similar to [H] we first perform a Mellin transform 

duo I a\z^ 



r{z,ai;p,a) = j 



a 



rM,P) (43) 



corresponding to Laplace transformation in rapidity. Here the ^-integration runs parallel to 
the imaginary axis to the right of the origin and, due to rapidity ordering, we use aiz"^ > a p^ . 
Eq. ()42|) becomes 

ujr^{zoi,p) = 6{p-Zoi) + ^ min<^^,l^ r^{z2i,p). (44) 

2vrp 2 Jo ^21 V Uoi Jy 



We perform another Mellin transformation 



dX fp"^^ 
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with the A-integration contour being the same as for the cj-integration above and assuming that 



z < r for all -Zj/s. We obtain 



which gives 



The solution of Eq. (ji^ is 



^r^x = — ^ + irT7 u^^^ 



1 A(cj-A) 
^-A = -3 ..r. .. X. ^.r (47) 



7rp2 tj[Acj- A2 - f] 



n J 2m 2m \ a p'^ J V^oi / uj\Xuj — X'^ — ^\ 

Similar to jH] we can first perform the integration over uj in Eq. ()48|): arguing that the 
high-energy asymptotics is dominated by the rightmost pole in uj plane and picking that pole 
yields 



TT p"^ J 2m \ a p^ J 2A^ + a 

Eq. (j49|) can be evaluated in the saddle point approximation near the saddle point at A = 
a/q;s/2, which gives 

r(^oi,«i;p,«) « ^ev^(^-^) f4) ' (50) 



^01 



leading to the well-known intercept of ^/2a~s |19| 1201 • A more careful evaluation of Eq. 
can be done by direct integration. First rewrite Eq. as 



r(2oi,ai;p, a) 



dX /ai\A ^ 1 



TT p'^ I 2m \ a J ^—^ n\ 

^ n=0 



«^ 1 ^01 
2A 



(51) 



2A2 + a 



The expression obtained in Eq. ()51|1 has poles at A = and at A = -^i^J'aJ2. From Eq. PH)) it 
is clear that the poles at A = -^i^fEi^j2 do not give leading high energy asymptotics. Concen- 
trating on poles at A = we neglect 2 A^ in the 2 A^ + Os denominator of Eq. (j3T|) . integrate each 
or the terms in the series by picking up the pole at A = and resum back the series obtaining 



r(zoi,ai;p, a) 



TT p2 



ai 2,2 



\ 2 In^ 



«sln(^) / / /._2 



Now we turn to studying the behavior of the nuclear modification factor under quantum 
evolution. First we note that the cross section in Eq. ()38p can be further simplified assuming 
transverse translational invariance of the solutions of the non-linear evolution equation for the 
dipole scattering amplitude, Eq. ()36|) . which can be achieved by considering an infinite and 
homogeneous nucleus in the transverse plane. Under this assumption the dipole scattering 
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amplitude becomes a function of a single spatial variable, N{r,b,y) N{r,y), and, shifting 
the integration variables in in the following way: 

x — z = v, (53) 

y — z = w, (54) 

z, = b, (55) 

Eq. (jHHj) can be rewritten as 



(27r)^ 2-11 ai J J 



|2 



d^kdy'-'"' (27r)2 27r ai 7 V 

[iV(|t;|, y) + N{\wly) - N{\v- w\,y)] . (56) 

Noticeably, the integrations over the arguments of the quark probability density r obeying 
the linear evolution equation ()32j) and the one over the arguments of the dipole scattering 
amplitude decouple. Thus, the linear evolution ()32j) only introduces to the cross section a 
rapidity-dependent prefactor, which carries no information at all of the nucleus and is therefore 
irrelevant for the nuclear modification factor. This prefactor can be calculated using the ap- 
proximate solution of the reggeon evolution equation given in Eq. (jH21)- The limits of integration 

2 

are set by the condition yi — y > In ^^^^ = yi — y > 0, yielding 



J d^zr{zjQ,z^,ai;z,a) = ^^^a^ h (^^2as{yi-y) ^ • (57) 

Finally, after performing analytically three of the integrals in the Fourier transform, Eq. (j56|) 
turns into 

dcr"" , . f ,2, yi-y T f L - , X , "1 



(P'k dy 



liO'"i) = y d% (^^2 as {yi -y)\n^^ 



dv 



2 1 
— Jiikv) — V Jo(k v) In — - 
k V A 



N{v,y). {51 



We now evaluate the nuclear modification factor RpA at increasing values of rapidity towards 
the proton fragmentation region, where the contribution to the cross section given in Eq. ()58p 
is the dominant one. The dipole scattering amplitude entering this equation is obtained by 
numerically solving Eq. (|36p. with initial conditions given by the classical limit, Eq. (0), with 
Q'^A = 1 GeV^ for the nucleus and Qlp = 0.1 GeV^ for the proton, and putting as = 0.2. 
The contribution to the cross section corresponding to the valence quark coming from the 
nucleus is included using the quasi-classical formula (PT|) at all values of rapidity. This is a 
good approximation for forward enough rapidities, where this contribution is vanishingly small 
anyway. The results for RpA at rapidities AY = y — Y/2 = 0,1,2,3, 4, are shown in Fig. [7| 

Our results show that the Cronin peak present in the initial condition is wiped out by 
the evolution, turning into a relative suppression at more forward rapidities. The rate of 
disappearance of the initial enhancement is similar to the one found in [30j for inclusive gluon 
production. This is not an unexpected result, since the effects of the reggeon evolution in the 
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Figure 7: Evolution of RpA_{kT, y) for valence quarks with rapidity. The initial condition corre- 
sponds to = 1 GeV^ for a nucleus and Q^p = = 0.1 GeV^ for a proton. The rapidity 
interval is measured from mid-rapidity, AY = y — Y/2. 

pp and pA cross sections are divided out in the nuclear modification factor, whose rapidity 
dependence is entirely driven by the dipole scattering amplitude, as can be seen from Eq. (j^Hl)- 
Moreover, the nuclear modification ratio at large rapidity increases with increasing transverse 
momentum, and eventually approaches unity at asymptotically large momentum. The slightly 
irregular behavior of the nuclear modification factor plotted in Fig. U\ at large momentum is 
due to the numerical uncertainty associated to the Fourier transform in Eq. (jKH|) and has by no 
means any physical origin. 

3.3 Valence Quarks from the Nucleus 

Now we want to include quantum evolution corrections into the expression for valence quark 
production cross section given in Eq. corresponding to the case (b) described above, in 

which the valence quark originates in the nuclear wave function. As can be seen from Fig. ^ 
the inclusion of gluon emissions with rapidities greater than the rapidity y of the produced 
valence quark can be straightforwardly accomplished using the linear dipole (BFKL) evolution, 
similar to how it was done in [25J for gluon production. 

We start by rewriting Eq. (pT|l for the case of dipole-nucleus scattering, similar to how it 
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was done above in Sect. 3.1 for the case (a): 



d^k dy 



1 2a, Cf 



-y+Y/2 



.01) 



(27r)2 TT 1-^ ^\k\ 



{x - yf 



-{lL-y?QL ln{l/|x-?;|A)/4 



d^xd^yd^Be'^^^-y) ^ (-1)^+-'' 
--{^-yfQ%Mi/\x-y\^)/i 



{x - bj) ■ {y - h.) 
\x — 6jp \y — 

(59) 



Here and hi are transverse coordinates of the quark and the antiquark in the incoming dipole 
and B_ = (1/2)(6q + b^) is the transverse coordinate of the center of the dipole. To include 
the quantum evolution corrections in the rapidity interval between y and Y we first denote the 
number of dipoles with transverse coordinates ^qm^i' cit rapidity y generated by the evolution 
from the original dipole 6o;^i having rapidity Y by ?^i(&0)^i;^o''^i'j ^ ~ v)- '^^^^ quantity is 
determined from the following evolution equation 



dni{bQ,bi;bQ,,bi,;y) a^N^ 



dy 



2 7r2 



d%2 



"01 
"20 "21 



^1 (^0, ^2 ; ^0' ,by]y) + ni{h2,bi-hj^,,hy-y) 



-ni{bQ,bi]hj^,,hy]y) 



with the initial condition 



ni{k^,bi]bQ,,by-y = Qi) = Sib^ - b^,) 5{bi - b^,). 
The inclusion of harder gluon emissions is accomplished by the substitution 



da'' 



d^k dy 



01 J 



d%Q> d%y 6o', &1'; Y - y) 



da'' 



d^k dy 



(^O'l'y 



(60) 



(61) 



(62) 



Inclusion of softer gluon emissions in the rapidity interval between and y is somewhat 
more involved. The structure of the quasi-classical result in Eq. (j?H) crucially depends on the 
fact that the t-channel quark exchanges in Fig. El are instantaneous. Attaching a single rung of 
Reggeon evolution to the quark lines in Fig. IHl would immediately separate the quark lines into 
a long-lived s-channel part and the short-lived t-channel part, as shown in Fig. |H1 Now the s- 
channel quark line allows for inclusion of quantum evolution corrections as was done in Eq. (jH^ 
and in Fig. IHl Thus it is the diagram in Fig. |H1 which gets dressed by quantum corrections. The 
answer for valence quark production cross section will consist of gluon production term given 
by the diagram in Fig. El plus the diagram in Fig. |H1 which should be enhanced by quantum 
corrections. Indeed, the emissions of gluons with rapidities larger than the rapidity y of the 
produced quark are allowed in Fig. |21 and do not significantly modify its structure. The real 
difference between Figs. ElandlHlis due to evolution between the produced quark and the target 
nucleus, which is allowed only in the latter. In fact, the quasi-classical term from Eq. ()21|) and 
Fig. El will become negligible compared to Fig. |H1 when quantum corrections are included. The 
situation is similar to the interactions of a current for dilaton-gluon coupling with the nucleus, 
as was first considered in j^. Since our expressions here should work for very high energies, 
we will neglect the contribution of Fig. El compared to that of Fig. El after quantum evolution 
is included. 
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Figure 8: Valence quark production in the case when the valence quark originates in the target 
nucleus with one rung of small-x evolution included. This diagram in turn serves as a lowest 
order diagram for the inclusion of small-x evolution corrections discussed in the text. 



The calculation of the diagram in Fig. |H1 proceeds along the lines of |60j. The wave function 
squared for an incoming quark at transverse coordinate decaying into a gluon which then 
decays into a hard massless quark and a soft anti-quark is given by jUOl 

/ a.\'^ X ■ (x — z) y ■ (y — z) + en Xi Zj % zi 
y;z,a) = ACJ-] ' , (63) 

The hard quark (to be produced) has transverse coordinates x and y to the left and to the 
right of the cut correspondingly, as shown in Fig. |H1 A much softer anti-quark has a transverse 
coordinate z_ and a fraction a -C 1 of the gluon's longitudinal momentum. Also, ei2 = 1, 
£21 = -1 and en = €22 = 0. 

Since the emissions of harder gluons which are described above in the dipole model frame- 
work lead to production of a dipole with transverse coordinates which for the moment we will 
denote by 6q and for the quark and for the anti-quark, we rewrite the wave function from 
Eq. (|63|) for the case of an incoming dipole (instead of incoming quark) as 

1 

^(^, m z,a) {-IY+^ ^{x -bi,y- by, z, a) (64) 

i,j=0 

where 

Hx-k,y-by,z,a)=AC,(^^^ (, _ y . _ ,^.)' (, _ _ 



X [(^ - hi) -{x- z) (y- bj) ■ iy-z) + emn {x - b^m {x - z)n Gki {v " bj)k {y - z)i] . (65) 

The interaction of the soft quark with transverse coordinate z_ with the nucleus is limited 
to exchange of two t-channel quarks, as shown in Fig. |H1 Indeed, since this quark has the 
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same transverse coordinate on both sides of the cut, the multiple rescattering corrections (with 
gluon exchanges like the ones shown in Figs. Q and 121) cancel via real- virtual cancellations. The 
interaction of the quark at z_ with the nucleus is then given by (see Appendix B of |9j) 

r!°!U,«) = ^^i^pT^(6)lnA (66) 

where s = s e"^"*"^ is the center-of-mass energy of the incoming gluon-nucleus system in Fig. |H1 
with s the center-of-mass energy of the whole collision. A is some infrared cutoff and we took 
into account the fact that nucleon contains valence quarks. 

The hard quark having transverse coordinates x and y in Fig. |H1 also can interact with the 
nucleus. Since the transverse coordinates of the hard quark are different on both sides of the 
cut [x_ and y\ real- virtual cancellations do not take place, and the interaction of this quark 
with the target brings in a Glauber- Mueller factor of 

Combining Eqs. (j64|) . and (j67|) and adjusting the overall normalization similar to 

how it was done in ,60^ we write the following expression for the quasi-classical valence quark 
production cross section in Fig. |H1 

'^^V(6^,J = -J-^ [ (fxd%d''z(fBe'^<^~y) [ daJ2^i^-k^y-bj'^^^^) 



i,j=0 

X r(2(^, a) e"(^-I^)^Q^.i-(VI^-^|A)/4. (68) 

Eq. (jUH|) makes inclusion of quantum corrections into the diagram of Fig. |H1 quite straight- 
forward. The emission of harder gluons is accomplished by the substitution in Eq. (jU^ as 
discussed above. All the possible emissions of softer gluons are illustrated in Fig. ^ A gluon 
may be emitted and absorbed by the harder quark, leading to nonlinear small-x evolution of 
the quark-anti-quark dipole x, y. Similar to |j2^ such evolution is included in the expression 
(jUHj) by substituting 

e-i--yfQ%Hi/k-yW/^ _ i-N{x,y,y) (69) 

where N{x, y, y) is given by the solution of Eq. (jHE))- In making the substitution of Eq. we 
are including gluon emissions both before and after the interaction with the target pKj . 

A gluon may not be emitted by the hard quark at x or ?/ and absorbed by the soft anti- 
quark at z_ (and vice versa) because such diagrams are Ai'c-suppressed. The only other allowed 
s-channel gluon lines are the emissions and absorptions of the gluon by the soft anti-quark 
at z_ as is also shown in Fig. IHl The gluon emitted this way can only be emitted before the 
interaction with the target: late-time emissions cancel by real- virtual cancellations [HS]- Since 
these gluons are emitted before the interaction on both sides of the cut, and as their transverse 
momenta are integrated over, their interactions with the target nucleus will cancel, just like 
they did for the anti-quark at z_. The result of such gluon emissions would thus be a linear 
evolution similar to the one shown in Fig. IHl but with rf^i{z,a) from Eq. as the initial 
condition. The evolution equation reads 



rint{x, y; z, a) = rl^t{z,a) + — / ^ r^U, z; Z2, a ) . (70) 



2 Jag « V± ~ ±2) 
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nucleus 



Figure 9: Some examples of quantum evolution corrections to the diagram in Fig. |H1 one may 
have emission and absorption by the hard and soft quark lines, as shown in the figure and 
discussed in the text. 

Just like for Eq. ()32|) the upper limit of a'-integral is derived by requiring ordering in rapidity, ao 
is some initial value of a at which the initial conditions are set. In the first step of the evolution 
in Eq. (|7n|l the ordering in rapidity is imposed by using both distances \z_ — x\ and \z — y\ as 
cutoffs, since the hard quark is the one we are tagging on and it has two different transverse 
coordinates on both sides of the cut. The subsequent steps of the evolution are easier, since 
there the harder quark will only have one coordinate. The resulting quantity rint{w_,w] z,a) 
obeys the following evolution equation, which is obtained from Eq. (fTUj) by simply putting 
x = y = w: 



(0)/ N , (^s 

Z TT 



-minU,^j da' d'z2 



a [Z-Z2) 



^2 ^int\Z_, Z_, Z_2, CM 



(71) 



With the help of Eqs. (jH^ . (fH^jl and (|7n|l we write the expression for the valence quark 
production cross section with the quark originating in the nuclear wave function including the 
effects of small-x evolution 



d^k dy 



(^oi) = 2 (2^7r)4 / d.^^ d^^o' d^bi' ni{bQ,b^;bf),,b^,;Y - y) e 



ik-{x—y) 



X 



nl 1 



ix-bi,,y- b.r,z, a) nmix, y; z, a) [I - N{x, y, y)]. 



(72) 



Integrating Eq. (f72j) over kx gives the following expression for the rapidity distribution of 
the stopped baryon number 



da 
dy 



2(2 



— -r / d^xd'^zd^Bd%Qid^biini{bQ,b^]bQ,,b^,]Y — y) 
Iny J 
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X da ^ix — bi,,x — bj,]z,a)rint{x,x;z,a). (73) 

It is interesting to note that, as could be expected, all the non-linear effects disappeared from 
da^/dy. 




a nucleus b 



Figure 10: Fan diagrams corresponding to the 'a' and 'b' contributions to the valence quark 
production cross section including the quantum evolution effects given by Eqs. (jHHj) and (f?^ 
corresp ondingly. 

The essential ingredients of Eq. (f?^ are the Reggeon exchange amphtude Tint obeying the 
linear equation (fTOj) summing up powers of as In^ s, and the amplitudes ni and obeying Eqs. 
(jUUI) and (jHUj) correspondingly, which sum powers of as In s. Note that at high kr, in the regime 
where evolution is linear, one can neglect in Eq. ()72j) . and, keeping only the linear (BFKL) 
part of Eq. for in Eq. (jHEj), one can show that each of the equations (PH|1 and (f?^ 
become dependent on a convolution of the BFKL pomeron and Reggeon exchange amplitudes. 

The fan diagram representation of the valence quark production cross sections calculated in 
Eqs. fj38|) and ()72j) are shown in Fig. ^1 The diagram 'a' in Fig. illustrates a contribution 
to Eq. ()38p with the valence quark originating in the proton: the diagram has a linear Reggeon 
evolution between the proton and the produced valence quark, with a full non-linear evolution 
of Eq. (|^ represented by fan diagrams between the produced quark and the target nucleus. 
The diagram 'b' in Fig. displays a contribution to Eq. (f?^ with the valence quark originating 
in the nucleus: it has a linear BFKL evolution between the proton and the produced valence 
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quark, a linear Reggeon evolution between the produced quark and the target, along with a 
possibility of a non-linear splitting attached to the produced valence quark line. 

4 Conclusions 

Above we presented calculations of valence quark production far from the fragmentation region 
in pA collisions. The answers for production cross section in the quasi-classical approximation 
are given by Eqs. ()2()|1 and ^IT^ corresponding to the valence quark originating in the proton 
and nucleus. The quasi-classical results lead to Cronin-like enhancement of nuclear modification 
factor for stopped baryons, as shown in Fig. Obtained transverse momentum spectra of the 
produced valence quarks scale as ~ at high kx and at fixed rapidity, indicating higher 

sensitivity of quark production to the physics at the ultraviolet end of the spectrum. 

Quantum evolution corrections, in the sense of double logarithmic reggeon evolution of 
[m I2DI 1211 1221 and single logarithmic gluon evolution of jTHl 1^ are included in the valence 
quark production cross section with the answers given by Eqs. and (f7^ corresponding to 
the valence quark originating in the proton and in the nucleus. Eq. fj38|) contains the effects 
of double-logarithmic linear reggeon evolution in the rapidity interval between the produced 
quark and the projectile along with non-linear gluon evolution effects between the quark and 
the target nucleus. The nuclear modification factor of stopped baryons coming from the pro- 
jectile proton (due to Eq. (jHHjl ) is explored and shown to give suppression of the produced 
baryons at higher energies and/or rapidities (see Fig.[7j). Therefore, based on the valence quark 
production spectrum and barring the effects of fragmentation functions and of non-perturbative 
mechanisms of baryon stopping [U El El 13 E] , we expect the stopped baryon spectrum to exhibit 
suppression in the forward rapidity region in d + Au collisions at RHIC. Indeed similar effects 
should be expected at LHC (possibly even at mid-rapidity WE\) if a pA run is to be conducted 
there. 

Eq. (j72p gives us the baryon stopping cross section with the valence quark coming from the 
nucleus. It includes linear BFKL evolution between the produced quark and the projectile and 
a combination of linear reggeon evolution and the non-linear evolution from Eq. (j36|) in the 
rapidity interval between the quark and the target. The study of the impact of small-x reggeon 
evolution in Eq. ()72p on nuclear modification factor is left for future work. Qualitatively we 
expect weaker suppression than for baryons coming from the proton, but the onset of this 
suppression should happen at smaller rapidities. Indeed for the non-linear evolution effects to 
become important the valence quark should be sufficiently far away in rapidity from the target. 
Due to the double-logarithmic nature of reggeon evolution one should expect it to manifest 
itself at smaller rapidities than the gluon evolution. At the same time, as the valence quark 
production cross section in Eq. ()72j) falls off approximately as ~ e~^, the non-linear evolution 
effects in it are likely to become important in the rapidity region where baryon stopping in the 
nucleus (i.e., the cross section itself) is small. This is likely to be the case at RHIC, where the 
small-x evolution effects do not become important until forward rapidity: there baryon stopping 
will be dominated by the valence quarks coming from the projectile deuteron (or proton) and 
baryon stopping due to Eq. will prevail over that from Eq. (f72|) . In other words evolution 
effects in baryon stopping in the nucleus for d + A collisions at RHIC may become important 
only in the rapidity region where baryon stopping is dominated by baryons from the deuteron. 
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Therefor, in order to test Eq. ()72|) one has to perform pA colhsion experiments at the LHC, 
where even at mid-rapidity one would expect significant small-x evolution effects [IHI- There 
one may find a window in rapidity (probably close to mid-rapidity but also not too distant 
from the nuclear fragmentation region) where baryon stopping is dominated by the mechanism 
given by Eq. (f72j) with the non-linear evolution effects being important. 

The results of our calculations are generally applicable to pA collisions to quantify the 
amount of purely perturbative baryon stopping. In the high-p^ sector, for pt > Qs, our results 
are also applicable to nucleus-nucleus collisions. It is likely, though by no means proven, that 
the obtained cross sections would provide a qualitatively correct description of the nuclear mod- 
ification factor for the perturbative initial-state baryon stopping in nucleus-nucleus collisions 
even at pr < Qs- Indeed at rapidity close to the fragmentation region of one of the nuclei in 
y4 + y4 collisions one may argue that the saturation scale of that nucleus {Qs2) would be small, 
much smaller than the saturation scale of the other nucleus [Qsi), making the particle pro- 
duction processes describable in the framework of pA collisions for a wide range of transverse 
momenta pt > Qs2- Therefore in that rapidity region baryon stopping can be described by the 
Eq. ()38|) above. The amount of baryon stopping not far from the fragmentation region is large, 
allowing for an easy experimental verification of our results. 

We conclude by reiterating that suppression of high-p^ net baryons in p{d) + A collisions at 
forward rapidity at RHIC, as predicted above, would provide an independent test of the physics 
of Color Glass Condensate. (Non-perturbative baryon stopping mechanisms are not likely to be 
important at high-pj'.) Similar experiments can be carried out at the LHC. The proposed test 
of CGC at RHIC with valence quarks may not be as clean as electromagnetic probes jHS], which 
are free from uncertainties introduced by fragmentation functions, but it has the advantage of 
being possible to perform by analyzing the current experimental data generated at RHIC. 
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Appendix A 

In this Appendix we derive the small-x quark and gluon wavefunction of a fast quark moving 
in the light cone "plus" direction. The corresponding diagram is shown in Fig. ^2 At lowest 
order in perturbation theory it is given by 



where p and a are the momentum and helicity of the incoming quark, k and a' those of 
outgoing quark and a and A the color and polarization of the emitted gluon. The transition 
matrix element is given by 



{f\Hj\t)/Vk^- 



,+ 



(Al) 



{f\Hi\€)/./^ = gV 



a 




Ua'{p) 



(A2) 
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where the gluon polarization vector in in the A'^ = gauge is* e{q) = (0, -p-,e'^) with e'^ = 
-j^{l,iX) and A = ±1. 




Figure 11: Diagram for soft quark wavefunction calculation. 
Defining a = /p~^, the energy denominator is 

We note that ^ 



p+«(l — a) 



1 



2b° 



ip-^ 



and 



6" 



(A3) 



(A4) 



(A5) 



/c+ p'' 

where e-^' the totally antisymmetric tensor in two dimensions. Using io^e^^e^ = aXk ■ one 
immediately gets 



ij':y,{p,k, «) = 9^1 + « - (1 - a)a\r-j^^-^^6--' 



(A6) 



Denoting by Xi the transverse coordinate of the incoming quark and by X2 and the transverse 
coordinates of the outgoing quark and gluon respectively and Fourier-transforming to coordinate 
space one gets 



^^.'Afei, ai23, «) = <[1 + a - (1 - a)aX]^6{x^, + ax2^)^^K^6"''' 



^23 



(A7) 



Taking the limit a — > and integrating over the transverse position of the hard gluon sets 
^3 = ^1 due to the delta-function and yields 



2tt xl^ 



5° 



(A8) 



To define the small-a quark distribution function of a nucleon we multiply the previous wave- 
function by its complex conjugate evaluated at different transverse position of the soft quark, 



*Here we follow the convention in which the product of two arbitrary vectors is given by a • 6 = i(a+& + 
a~b'^) — a ■ b 
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average over the helicity of the initial quark and sum over the polarizations and colors of the 
gluon and integrate over the position of the initial quark (xi) obtaining 

dn'^ 1 a s-^ a (x — z) ■ (y — z) 

(T,A,a — 

as given in Eq. 0. The factor a/2 arises from the integration over the phase space of the 
gluon: da/2{l — a) ^ ady/2. 

Taking the limit 1 — a ^ in Eq. ()A7jl and proceeding analogously one derives the soft gluon 
distribution function of a valence quark 

Kyx{xu,c^) = 2ge^^^5^^' (AlO) 
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and 



dn^ 1 Iv^/.a / ^ r ^^ a., {x - z) ■ {y - z) 

d^zdy An 2 ^ ^' -n {x - zy{y - zf ^ ' 



a.X.a 



Appendix B 

In this appendix we derive the Gq qG scattering cross section in the high energy limit. The 
calculation is performed in the A'^ = gauge, with gluon polarization vectors taken in the 
same gauge e"*" = 0. The amplitude for the process is shown in Fig. [Hand is equal to 

M = uAk + l){igt^Y)j^{wtVW{p)4(k)ef{p - I). (Bl) 

r + ze ^ 

Squaring the amplitude, averaging over the helicity of the incoming quark and the polarization 
of the incoming gluon and summing for all possible helicities, polarizations and colors in the 
final state, one gets 



(iMp) = 1 — - — iMp 



^ ' a,b,cr,CT' ,X,X' 

= (2^)'^ (/2 f,,)2 t^[(^ + ^)rhy7'md,a{k)d^p{p - I) (B2) 

where d^^ik) = J2x^^(k)e'^{k) = ~g^a + is the gluon polarization sum in the light 

cone gauge. 

Keeping only the leading terms, i.e. those proportional to the center mass energy of the collision, 
s = 2p~k~^, and putting Ri — we get 

{\M\^)=4{2nr^^. (B3) 

Finally, dividing by the flux factor and integrating over the phase space of the produced quark 
and gluon one gets the final result 

G,->gG _ I ffi .g^^ 
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V 

Figure 12: Gq^qG scattering amplitude. 
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